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Abstract. We develop relative oscillation theory for Jacobi matrices which, 
rather than counting the number of eigenvalues of one single matrix, counts the 
difference between the number of eigenvalues of two different matrices. This is 
done by replacing nodes of solutions associated with one matrix by weighted 
nodes of Wronskians of solutions of two different matrices. 



1. Introduction 

Oscillation theory for second-order differential and difference equations has a long 
tradition originating in the seminal work of Sturm from 1836 [5]. Since then the 
subject is continuously growing and many monographs have been devoted entirely 
to this subject. The most recent one being the monumental treatise by Agarwal, 
Bohner, Grace, and 'Regan [Tj . One of the key results of classical oscillation theory 
is the fact, the A;'th eigenfunction has precisely k—1 nodes (i.e., sign flips) and for a 
suitably chosen solution of the underlying difference equation, the number of nodes 
of this solutions equals the number of eigenvalues below a given value. Our aim is 
add a new wrinkle to this classical result by showing that the number of weighted 
nodes of the Wronskian (also known as Casoratian) of two suitable solutions of two 
different Jacobi difference equations can be used to count the difference between 
the number of eigenvalues of the two associated Jacobi matrices. 

That Wronskians are related to oscillation theory is indicated by an old paper 
of Leighton [7] , who noted that if two solutions have a non- vanishing Wronskian, 
then their zeros must intertwine each other. However, it seems their real power was 
realized only later by Gesztesy, Simon, and Teschl in [3] with the corresponding 
extension to Jacobi operators given by Teschl [10] . For a pedagogical discussion 
we refer to the survey by Simon [8]. That these results are just the tip of the 
iceberg was discovered only recently by Kriiger and Teschl [4], [5], [6]. Our result 
generalizes the main result for the case of Sturm-Liouville operators from [5] to the 
case of Jacobi matrices. 
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To set the stage, let us fix some real numbers a(j) < 0, b(j), j — 1, • • • , N — 1 
and consider the Jacobi matrix 

\ 




b(N-2) a(N-2) 
a(N-2) b(N-l)J 

in the Hilbert space C^ -1 . Furthermore, let s±(z,n) be the solutions of the un- 
derlying difference equation (set a(0) = a(N — 1) = a(N) = —1, b(N) = 0) 

(1.2) a{n)u{n + 1) + b(n)u(n) + a(n — l)u(n — 1) = zu(n), n = 1, . . . ,N, 
corresponding to the initial conditions 

(1.3) s-(z,0)=0,s-(z,l) = l, s + (z,N) = 0,s + (z,N + l) = l. 

Note that s_(A, n) (resp. s_|_(A, n)) will be an eigenvector of H corresponding to 
the eigenvalue A e 1 if and only if s_(A, iV) = (resp. s + (A, 0) = 0). We will 
abbreviate s(z,n) — S—(z,n). 

We call n a node of a solution u of (|1.2|) if cither 

(1.4) u(n) = or u(n)u(n + 1) < 0. 
We say that a node no of u Zies between m and n if either 

(1.5) m < hq < n or no — m but it(m) ^ 0. 

#(m,n)( u ) denotes the number of nodes of u between m and n and #(m) = #(o,iv)(' u )- 
Then we have the following classical result alluded to before (see e.g., [5], [TTj): 

Theorem 1.1. Let H be a Jacobi matrix and s(z,n) a corresponding solution 
of the underlying difference eguation (|1.2|) corresponding to the initial condition 
s(z, 0) = 0. Then for every Asl the number of nodes of s(A, n) equals the number 
of eigenvalues of H below X: 

(1.6) #(,s(A))=#{£e,(i/)|E<A}. 

Here o~(H) denotes the spectrum of H , that is, the set of eigenvalues. 

To generalize this result we will now consider two Jacobi matrices Hq and H\ 
associated with the coefficients ao(n) = a>\(n) = a(n) and bo(n) respectively b\{n). 
The corresponding solutions will be denoted by Sj t ±(n), j — 0,1, in obvious no- 
tation. Given two solutions Uj of the difference equations associated with Hj we 
denote by 

(1.7) W n (uo,ui) = a(n)(u (n)ui(n + 1) - u (n + l)ui(n)) 

their Wronskian. As already anticipated we will relate the number of nodes of such 
Wronskians to the difference between the eigenvalues of Hq and Hi. Since this 
difference is a signed quantity, we will need to weight the nodes according to the 



/6(1) o(l) 

o(l) b(2) 


a(N-l) 

V o o o 
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sign of the difference between Hq and H\ as follows: Set 

if b {n + 1) - bi{n + 1) > and 
1, either W n (uo,Ui)W n+ i(uo,Ui) < 

or W n {u ,Ui) = and W n+ i(u ,Ui) ^ 0, 

(1.8) #n{u a ,u 1 ) = { if b Q (n + 1) - b 1 (n + 1) < and 

-1, either W n (u , ui)W n+ i(u , m) < 

or W n (u ,ui) ^ and W n+ i(u ,ui) = 0, 

t 0, otherwise. 

Then we say the Wronskian has a weighted node at n if f£ n (u ,ui) ^ 0. The 
number of weighted nodes of the Wronskian between and N is denoted as 



N- 



' ' - if II o ( 'in- a i. ! -f 0. 

With this notation our main result reads 



(i.9) #k,«i) = ]r#,K,ui)-<T' 



if Wo(-u ,mi) = 0. 



Theorem 1.2. Let iJo, -Hi be two Jacobi matrices with ao — a\ and Sj^±{z,n), 
j = 0,1, the corresponding solutions of the underlying difference equations. Then 
for Ao,Ai £ M the number of weighted nodes of W(sq^(Xq), si i +(Ai)) equals the 
number of eigenvalues of Hi below Ai minus the number of eigenvalues of Hq below 
or equal to A : 

#( S0) -(A ), si,+ (Ai)) = #(so,+(A ), si,-(Ar)) = 
(1.10) = #{S e a(Hx)\E < A x } - #{£ e a(H )\E < A }. 

Here o~(H) denotes the spectrum of H, that is, the set of eigenvalues. 

The proof is based on Prufer angles to be investigated in Section [2j It will be 
given in Section [3] 

An extension to Jacobi operators on N respectively Z is in preparation. 

2. Prufer angles 

Since any nontrivial solution of (|1.2| cannot vanish at two consecutive points 
we can introduce Prufer variables (p u (n),9 u (n)) in the usual way (cf., e.g., [Ill 
Chap. 4]) via 

(2.1) u(n) = p u (n) sin(0 u (n)), u(n + 1) = p u {n) cos(0 u (n)). 

Note that p u (n) > for all n 6 Z and 9 U (n) is only defined up to an additive integer 
multiple of 2ir, depending on n. For our further investigations it is essential to gain 
unique values for the Prufer angle and therefore we fix 9 U (0) and require 

(2.2) rWAl < \0u{n + 1)/tt] < [e u (n)/7r] + 1, 

where [x~\ = min{n 6 Z | n > x} denotes the usual ceiling function. Then the 
following easy result is well-known. 

Lemma 2.1. Define k, 7, T via 

(2.3) 9 u (n) = fc7r + 7, 6 u (n + 1) = kn + T, 7 e (0, tt], T e (0, 2tt], k e Z. 
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Then 

(0; f] iff n * s n °t a node, 



(2.4) 7 G 
and 

(2.5) r g 

Moreover, 



I (f , 7r] i/f n is a node, 

(0, 7r] jj(f n is not a node, 
(it, 2tt) iff n is a node. 



(2.6) 9 u (n) = kir+ | O 0„(n + 1) = (k + 1)tt. 

As a consequence we obtain 
Corollary 2.2. We have 

^ ^ ^ 9 u (n+ 1) _ J r^^l +1 if n is a node, 

7T 



|- flu (w) i otherwise. 



In particular, we obtain 



( , 8) #M . ^ _ L M0)j _ L 

7T 7T 

where = max{n G Z | n < a;} is the usual floor function. 

To find the analogous formula for the number of weighted nodes of a Wronskian 
we observe 

(2.9) W n {u ,ui) = -a(n)p Uo (n)p Ul (n)sm(A Uo , Ul (n)), 
where 

(2.10) A uo . ul (n) = 6 ui (n) - 9 U0 {n). 
Furthermore, note 

(2.11) W„+i(u ,wi) - W n (uo,ui) = (6 (n + 1) - b x {n + l))u {n + l)u x (n + 1). 
As a straightforward consequence of Lemma 12.11 we obtain 

Lemma 2.3. Fix some n and let 6j (n) = kjir + with 7j G (0, it] and 9j (n + 1) = 
kjir + Tj with Tj G (0, 2n] for j = 0, 1. TTien we /icwe 

(2.12) A UOjUl (n) = (fci-fc )7r+7i-7o cmd A„ 0iUl (n+1) = (fci-fc )7r+ri-r , 
w/iere 

(1) : either u and U\ have a node at n or both do not have a node at n, then 

(2.13) 7i-7oG(~,|) and T x - T G (-it, tt). 

(2) : ?ii ^as node a£ n ; but uq has a node at n, then 

(2.14) 7i - 70 G {-n, 0) and T x - T G (-2tt, 0). 

(3) : ui /ifls a node at n, &Tii uo has no node at n, then 

(2.15) 7i-7oG(0,tt) and T x - T G (0, 2tt). 
Now we are able to show 
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Lemma 2.4. Fix some n. Then, if bo(n + 1) > b\{n + 1), we have 

(2.16) rA U0 , Ul (n)M < \A UOiUl (n + l)/if\ < \A UOiUl (n)/*] + 1 
and if bo(n + 1) < b\{n + 1), we have 

(2.17) [A Uo , Ul (n)/7rl - 1 < \\ l0 , Ul (n + 1)/tt"| < rA Uo , Ul (n)/7r| . 



Proof. We will use the notation from Lemma 12.31 where we assume fco = fci = 
without loss of generality. In particular, Lemma 12.31 implies 



«o,«i ( n )/ 7r l ~ 1 — [A U0!U1 (n + 1)/tt] < |~A Uo>Ul (n)/7r] + 1. 

Hence, to show (|2.16p there are two cases to exclude. Namely, (i) A„ 0iUl (n) <E (0, ■§), 
A UOlU1 (ra+l) G (-tt,0] (from case (1)) and (ii) A Uo , Ul (ra) E (-7r,0), A M0;Ml ( rt + l) S 
(— 2ir, — 7r] (from case (2)). But in case (i) we obtain a contradiction from (|2.11|) : 

Wn + i(u ,tti) = Wn(no,m) + (b (n + 1) - bi(n + 1)) it (n + l)m(n + 1) . 
<o >o >o >o 

Similarly, in case (ii) equation (|2.11[) implies 

Wn+xjuo,^)^ = W n {uQ,Ui) + (b (n + 1) - bi{n + lWn + l)zti(n + 1) . 

>0 <0 >0 <0 

Equation (|2.17p can be established in a similar fashion. □ 

Lemma 2.5. Let n 6 Z, t/ien 

(1) : W^ n (tt ,«i) = W n+ i(« ,ui) =0 orWn^o^ijWn+^Mo.Ui) > implies 
/ 2 lg s p A n0ittl (n+ 1) ^ = ^ A M0 , Ml (n) ^ 

7T 7T 

(2) : Wn(uo,Ui)W n 4-i(uo,«i) < implies 

A U0 , Ul (n + l) = /r A " 0, ; i(n) l + 1 - «/6o(n + l)>6i(n+l), 
U ^gfflW ]-!, if6 (» + l)<6 1 (n+l). 



(2.i9) r 



(3): W n (uo,Ui) — and W n +i(uo,ui) ^ implies 

(2 20) r Auo, Ul (n + l) _ [r ^^ l+l, i/6o(n + l)>6i(n+l), 

U A "°: i(n) l, ifb (n + l)<h(n+l). 



(2.2i) r- 



(4): W n (u( h ui) cmd W n +i(uo,Ui) = implies 

.A UOtU1 (n + l) 1 _ f r A "°-; i(n) l, ifb (n + l)>b 1 (n+l), 

,' Au °: i( "' l -1, i/6 (n + l)<6i(n + l). 

A/oie i/iai in i/ie cases (2) -(4) we necessarily have bo{n + 1) =/= bi(n + I) . 

Proof. We will use the notation from Lemma 12.31 where we assume fco = fci = 
without loss of generality. Moreover, interchanging uq and u\ using A Ul)tl0 = 
-A„ 0itll (n) and 

— \x\ if x £ Z, 

— \x] + 1 otherwise, 

we see that it suffices to show one case 6o(n+ 1) > b\(n+l) or bo(n+ 1) < &i(n + l). 
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Suppose W n {uo,ui) = W n +i(uQ 7 ui) = and W n (uo,ui)W n +i(uo,ui) > do 
not hold, then by (|2.11j) we have 

W n+ i{u ,ui) - W n (u ,ui) = {b (n + 1) - bi(n + l))u (n+ l)tii(n+ 1)^0 

and hence bo(n + 1) ^ b\{n + 1). 

(1) and (2). Suppose W n (uo,ui) = W n +i {uo, m) = 0, then by (|2.9|) we infer 

sin(A Uo , Ul (n)) = sin(7i - 70) = 0, sm(A U0)Ul [n + 1)) = sinfTi - T ) = 0, 

where 70, 71 G (0,7r]. Thus 70 = 71 and we have case (1) of Lemma 12.31 which 
implies Ti — To G (— 7r,7r) and we conclude T\ — To = 0. In summary, A Uo , Ul (n) = 
A UOiUl (n + 1) = as claimed. 

Next suppose W n (uo,Ui)W n +i(uo,ui) ^ 0, then by (|2.9|) the sign of the Wron- 
skian at n equals the sign of sin(A UOjUl (n)) and hence (|2. 16[) respectively (|2.17|) 
finish the proof of case (1) and (2). 

(3) . By (|2.9j) we conclude A U0)ttl (n) — 71 — 70 = mod 71", where 70, 71 € (0, 7r] 
and thus 71 — 70 = 0. So we have case (1) of Lemma [231 and hence A U0)ttl (n + 1) = 
r! -T e (-7T,7r). That is, 

rA UOiUl (n)/7rl < [A U0>Ul (n + l)/7r] < |"A„ , Ul (n)/7rl + 1 

and (I2.17|) finishes the proof of case (3) for bo(n + 1) < b\{n + 1). 

(4) . By (|2.9|) we have A UOjUl (n + l) = Ti — To =0 mod it and Lemma l2~3l leaves 
us with the following possibilities 

(a) A UOiUl (n) e (-|,f) and A Uo , Ul (n + 1) = 0, 

(b) A„ 0iUl (n) e (-7T,0) and A„ 0)Ml (n + 1) = -7T, 

(c) A tt0iUl (n) G (0,7r) and A„ 0iUl (n + 1) = ir. 

and (j2~T6)) shows (4) if b Q (n + 1) > 61 (n + 1). □ 

As a consequence we obtain the desired formula 
(2.22) #(uq, Ul ) = rA„ 0)U1 (iV)/7rl - [A U0<U1 (0)/ttJ - 1. 

3. Proof of the main theorem 

Our strategy will be to interpolate between H and H\ using H e = (1 — e)i?o + 
£-ffi, that is, a e (n) = a(n) and 6 e (n) = (1 — e)bo(n) + ebi(n). If u e is a solution of 
the difference equation corresponding to H e , then the corresponding Priifer angles 
satisfy 

3.1) e (n) = rT~27T ' 

a{n)pi(n) 

where the dot denotes a derivative with respect to e. 
Lemma 3.1. We have 



~ E m , n+ i(M m ) - h{m))s e ,+ {z, m) 2 , 
ELi(M ra ) - & i( m )K ,-{z,m) 2 . 



(3.2) W n (s e ,±(«),A B ,±(z)) = < 
Proof. Summing (|2.11[) we obtain 

LmdWl" 1 ) - b 1 [m))s £ _{z,m)Si -{z,m). 
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Now use this to evaluate the limit 

s s ,±(z) - s i,±( z ) 



\imW 7 



□ 



Denoting the Prater angles of s e .±(A, n) by e ,±(A, n), this result implies for 
b Q - bx > 0, 

a n \ EL+i(M m )-i'iWK,+(^'") 2 , n 

B ^ n) = (nK + (A,n)» " °' 

(3.3) 4 _(A,n) = - ^-^"W -^Wk-t^) 2 > Q 

«WPe ,-(A,nr 

Furthermore, we have the following result from classical perturbation theory. We 
add a simple direct proof for convenience of the reader. 

Lemma 3.2. Suppose bo — b\ > (resp. bo~b\ < 0). Then the eigenvalues of H e 
are analytic functions with respect to e and they are decreasing (resp. increasing). 

Proof. First of all the Priifer angles 9 £j ±(X, n) are analytic with respect to e since 
s e,±(A, n) is a polynomial with respect to e. Moreover, A 6 cr(H s ) is equivalent 
to e ,+ (A, 0) = mod 7T (resp. 9 e ^(X,N) = mod n) and monotonicity follows 
from' JO). □ 



In particular, this implies that P(H £ ) = f/={E G a(H £ )\E < A} is continuous 
from below (resp. above) in e if bo — b\ > (resp. bo — b\ < 0). 
Now we are ready for the 



Proof of Theorem M.SX It suffices to prove the result for #(so,+ (Ao), si,_(Ai)), where 
we can assume Ao = Ai = without restriction and set s St ±(n) — s e ,±(0,n) for 
notational convenience. We split 6o — b\ according to 

bo — b\ = & + — 6_ , b + , 6_ > 0, 

and introduce the operator H— = Ho — 6- . Then H _ is a negative perturbation of 
Ho and Hi is a positive perturbation of H- . 
Furthermore, define H s by 



H, 



Ho + 2e(H_-H Q ), eG [0,1/2], 

< H- + 2(e-l/2)(Hi-H-), £G [1/2,1]. 

Let us look at (using (|2~22| ) 

Q{e) = #( So ,+, fl e ,_) - [A £ (A0A1 - LA £ (0)/7rJ - 1, A e (n) = A S0;+ , Sc ,_ (n) 

and consider £ G [0, 1/2]. At the left boundary A e (0) remains constant whereas at 
the right boundary A e (JV) is increasing by (13. 3|) . Moreover, it hits a multiple of 7r 
whenever G a(H e ). So Q(e) is a piecewise constant function which is continuous 
from below and jumps by one whenever G o~(H £ ). By Lemma l3.2l the same is true 
for 

P(e) = #{E e o-(H e )\E < 0} - #{E G a(H )\E < 0} 

and since we have Q(0) = P(0), we conclude Q(e) = P(e) for all £ G [0, 1/2]. To 
see the remaining case £ = [1/2,1], simply replace increasing by decreasing and 
continuous from below by continuous from above. □ 
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